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Thermal fluctuations of quasiparticle number are included making use of the 
secondary Bogolyubov's transformation, which turns quasiparticles operators 
into modificd-quasiparticle ones. This restores the unitarity relation for the gen- 
eralized single-particle density operator, which is violated within the Hartree- 
Fock-Bogolyubov (HFB) theory at finite temperature. The resulting theory is 
called the modified HFB (MHFB) theory, whose limit of a constant pairing in- 
teraction yields the modified BCS (MBCS) theory. Within the MBCS theory, 
the pairing gap never collapses at finite temperature T as it does within the 
BCS theory, but decreases monotonously with increasing T. It is demonstrated 
that this non-vanishing thermal pairing is the reason why the width of the 
giant dipole resonance (GDR) does not increase with T up to T ~ 1 MeV. At 
higher T, when the thermal pairing is small, the GDR, width starts to increase 
with T. The calculations within the phonon-damping model yield the results 
in good agreement with the most recent experimental systematic for the GDR 
width as a function of T. A similar effect, which causes a small GDR width at 
low T, is also seen after thermal pairing is included in the thermal fluctuation 
model. 



1. Introduction 

It is well known that infinite systems undergo a sharp phase transition 
from the superfluid phase to the normal-fluid one at finite temperature T. 
Marked by a collapse of the pairing correlations (pairing gap), and a near 
divergence of the heat capacity at a critical temperature T c , this phase 
transition is a second-order one. The critical temperature is found to be 
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T c ~ 0.567A(0) for infinite systems, where A(0) is the pairing gap at zero 
temperature T = 0. 1 

The application of the BCS theory and its generalization, the Hartree- 
Fock-Bogolyubov (HFB) theory, to finite Fermi systems paved the way to 
study the superfluid-normal (SN) phase transition in nuclei at finite tem- 
perature. 2-4 Soon it has been realized that the BCS and HFB theories 
ignore a number of quantal and thermodynamic fluctuations, which be- 
come large in small systems because of their finiteness. As a consequence, 
the unitarity relation for the generalized particle-density matrix R, which 
requires R 2 = R, is violated. In deed, within the HFB theory at T ^ 0, one 
has Tr[R 2 {T) - R(T)] = 25N 2 = 2£.<SA/? - 2£\n;(l - m) > where 
rii = [c /3Ei + 1] _1 is the occupation number of non-interacting quasiparticlcs 
with energy Ei at temperature T = 1//3 on the i-ih orbital. Large ther- 
mal fluctuations smooth out the sharp second-order SN phase transition. 
As the result the pairing gap does not collapse as has been predicted by 
the BCS theory, but decreases monotonously as the temperature increases, 
and remains finite even at rather high T. 5 ~ 7 So far these fluctuations were 
taken into account based on the macroscopic Landau theory of phase tran- 
sitions. 5 ' 6 This concept is close to that of the static-path approximation, 
which treats thermal fluctuations on all possible static paths around the 
mean field. 7 

It will be shown in the first part of this lecture that the recently proposed 
modified-BCS (MBCS) theory, 8 ' 9 and its generalization, the modified-HFB 
(MHFB) theory 10 take into account the fluctuations of quasiparticle num- 
ber in a microscopic way. The MHFB theory restores the unitarity relation 
by explicitly including the quasiparticle-number fluctuations, making use 
of a secondary Bogolyubov transformation from quasiparticle operators to 
modified quasiparticle ones. In the limiting case of a constant pairing in- 
teraction G the MHFB equation is reduced to the MBCS one. 

The second part of the lecture represents an application of the MBCS 
theory in the study of the damping of giant dipole resonances (GDR) in hot 
nuclei, which are formed at high excitation energies E* in heavy-ion fusion 
reactions or in the inelastic scattering of light particles (nuclei) on heavy 
targets. The 7-decay spectra of these compound nuclei show the existence 
of the GDR, whose peak's energy depends weakly on the excitation energy 
E* . The dependence of the GDR on the temperature T has been experi- 
mentally extracted when the angular momentum of the compound nucleus 
is low, as in the case of the light-particle scattering experiments, or when it 
can be separated out from the excitation energy E* . These measurements 
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have showed that the GDR width remains almost constant at T < 1 MeV, 
but sharply increases with T up to T ~ 2 - 3 MeV, and saturates at higher 
T. 11 The phonon-damping model (PDM), proposed by the lecturer in col- 
laboration with Arima, 12 explains the GDR width's increase and saturation 
by coupling the GDR to non-collective particle-particle (pp) and hole-hole 
(hh) configurations, which appear due to the deformation of the Fermi sur- 
face at T 7^ 0. It will be shown that, by including non-vanishing MBCS 
thermal pairing, the PDM is also able to predict the GDR width at low T. 

2. Modified HFB theory at finite temperature and its limit, 
modified BCS theory 

2.1. HFB theory 

The HFB theory is based on the self-consistent Hartree-Fock (HF) Hamil- 
tonian with two-body interaction 

h = Y1 Ti o a \ a j + \ Yl v m a i a ] aiak ' w 

ij ijkl 

where .. denote the quantum numbers characterizing the single-particle 
orbitals, %j are the kinetic energies, and Vijki are antisymmctrized ma- 
trix elements of the two-body interaction. The HFB theory approximates 
Hamiltonian (1) by an indcpendent-quasiparticle Hamiltonian _ffjjFB 

H - flN w #hfb =E Q +J2 E M a i ' ( 2 ) 

i 

where N is the particle-number operator, \x is the chemical potential, Eq 
is the energy of the ground-state |0), which is defined as the vacuum of 
quasiparticles: 

ai\0) = , (3) 

and Ei are quasiparticle energies. The quasiparticle creation a\ and de- 
struction Qfj operators are obtained from the single-particle operators a\ 
and a>i by the Bogolyubov transformation, whose matrix form is 

with the properties 

UU^ + VV^ = 1 , UV T + VU T = , (5) 
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where 1 is the unit matrix, and the superscript T denotes the transposing 
operation. The quasiparticle energies Ei and matrices U and V are deter- 
mined as the solutions of the HFB equations, which are usually derived by 
applying either the variational principle of Ritz or the Wick's theorem. 

At finite temperature T the condition for a system to be in thermal 
equilibrium requires the minimum of its grand potential Q 

n = £ -TS - fiN , (6) 

with the total energy £ , the entropy S, and particle number N, namely 

5fl = . (7) 

This variation defines the density operator V with the trace equal to 1 

TrX> = l, Sfl/6V = (8) 

in the form 

V = Z- 1 e-^ H ~^\ Z = Tr[e-^-^)] , (i = T _1 , (9) 

where Z is the grand partition function. The expectation value -< O y of 
any operator O is then given as the average in the grand canonical ensemble 

< 6 y= Tr(D0) . (10) 

This defines the total energy £, entropy 5, and particle number N as 

£ = Tr(VH) , S = -Tr(X>hiD) , N = Tr(VN) . (11) 

The FT-HFB theory replaces the unknown exact density operator V in Eq. 
(9) with the approximated one, 2?hfb, which is found in Ref. 3 by substi- 
tuting Eq. (2) in to Eq. (9) as 

2?hfb = Hl^M + (1 - rn)(l-tfi)] , (12) 

i 

where Mi is the operator of quasiparticle number on the i-th orbital 

M = a\ai , (13) 

and rii is the quasiparticle occupation number. Within the FT-HFB theory 
rii is defined according to Eq. (10) as 

rn = (Afi) = , (14) 

where the symbol (. . .) denotes the average similar to (10), but in which 
the approximated density operator 2?hfb (12) replaces the exact one, i.e. 

(6) = Tr(2?HFBO) • (15) 
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The generalized particle-density matrix R is related to the generalized 
quasiparticle-density matrix Q through the Bogolyubov transformation (4) 
as 



R=U f QU , (16) 



where 



with 

U =( U V V u) > Utf = \. (18) 

The matrix elements of the single-particle matrix p and particle pairing 
tensor r within the FT-HFB approximation are evaluated as 

Pij = (ata,-) , T y - = (ajdi) , (19) 

while those of the quasiparticlc matrix q are given in terms of the quasi- 
particle occupation number since 

qij = (a]ai) = Sijrii , t y - = (aj-aj) = , (20) 

which follow from the HFB approximation (2). Using the inverse transfor- 
mation of (4), the particle densities are obtained as 3 

p = U T nU* + V\l-n)V , r = U T nV* + V\l - n)U . (21) 

By minimizing the grand potential according to Eq. (7), the FT-HFB 
equations were derived in the following form 3 

("MS)- 

where 

H = T-p + T, Tij =^2v ik:j ipi k , Ajj = ^y^VjjkiTki ■ (23) 

kl kl 

The total energy £, entropy S, and particle number N from Eq. (11) are 
now given within the FT-HFB theory as 

£ = Tr[(T +\T)p+\ Art], (24) 
S = - ^[ndnn; + (1 - n,)ln(l - m)\ , (25) 
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N — Trp , (26) 

from which one can easily calculate the grand potential Q (6). 

At zero temperature (T = 0) the quasiparticle occupation number van- 
ishes: rii =0, and the average (15) reduces to the average in the quasiparticle 
vacuum (3). The quasiparticle-density matrix Q (17) becomes 

Q(T = 0) = Qo = , for which Q 2 = Qo. (27) 

Therefore, for the generalized particle-density matrix i?o = i2(T = 0) the 
following unitarity relation holds 

Rl = Ro , where R = U ] Q Q U . (28) 

However, the idempotent (28) no longer holds at T ^ 0. Indeed, from Eqs. 
(16) and (17) it follows that 

R- R 2 =U\Q -Q 2 )U , (29) 

which leads to 

Tr(R-R 2 ) = Tr{Q-Q 2 ) = 2j2n t (l-n l ) = 2(5N) 2 ^0 , (T + 0) . 

(30) 

The quantity JA" 2 = — n i) m Eq. (30) is nothing but the 

quasiparticle-number fluctuation since 

5N 2 = (A -2 ) - (A/ - ) 2 = <E#i + E^> -E n '-E 

= E«i(i-«i) = E*- A/ ?' ( 31 ) 

where = n,(l — rij) is the fluctuation of quasiparticle number on the 
i-th orbital. Therefore, in order to restore the idempotent of type (28) at 
T / a new approximation should be found such that it includes the 
quasiparticle-number fluctuation in the quasiparticle-density matrix. 

2.2. MHFB theory 

Let us consider, instead of the FT-HFB density operator X>hfb (12), an 
improved approximation, f>, to the density operator V. This approximated 
density operator V should satisfy two following requirements: 
(i) The average 



((6))=Tr(Vd), 



(32) 
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in which V is used in place of V (or 2?hfb), yields 

R = U^QU 



(33) 



for the Bogolyubov transformation U (18), where one has the modified 
matrices 



R 



P r 
-f* 1 - p* 



Q 



with 



Pij = (( a } a *» ' 



q t 

-t* 1 - q* 



{{a jai )) , 



Qij = (( a ] a i)) = S i3 n i . tij = (( a J a i)) = A ij 



(34) 



(35) 
(36) 



instead of matrices R and Q in Eqs. (17), (19), and (20). The non-zero 
values of iij in Eq. (36) are caused by the quasiparticle correlations in the 
thermal equilibrium, which are now included in the average ((•••)) using 
the density operator f>. 

(ii) The modified quasiparticle-density matrix Q satisfies the unitarity 
relation 



(0) 2 = Q ■ 



(37) 



The solution of Eq. (37) immediately yields the matrix A in the canonical 
form 



A = \/n(l — n) 



(1 -Ai\ 
Ai 



-A 2 
A 2 



V 



A 4 = y/ni{\ ~ Hi) 



J 



(38) 

Comparing this result with Eq. (31), it is clear that tensor A consists of 
the quasiparticle-number fluctuation Sfifi = ^/fij(l — n»). From Eq. (33) it 
is easy to see that the unitarity relation holds for the modified generalized 
single-particle density matrix R since R — R 2 — W (Q — Q 2 )U = due to 
Eq. (37) and the unitary matrix U. 

Let us define the modified-quasiparticle operators a\ and on , which be- 
have in the average (32) exactly as the usual quasiparticle operators a] and 
ai do in the quasiparticle ground state, namely 



((a\a k )) = ((al4)) = ((a k ai))=0. 



(39) 



February 9, 2008 0:45 



WSPC - Proceedings Trim Size: 9in x 6in DangPredeal06 



8 



In the same way as for the usual Bogolyubov transformation (4), we search 
for a transformation between these modified-quasiparticle operators (at, 
cti) and the usual quasiparticle ones (at, a,) in the following form 

(40) 



w z 
z* w* 



with the unitary property similar to Eq. (5) for U and V matrices : w f + 
zzt = 1 . Using the inverse transformation of (40) and the requirement 
(39), we obtain 



««k)> 



k 



(41) 



From this equation and the unitarity condition (28), it follows that zz 



and 



,t 



1 — n. Since 1 — n and n are real diagonal matrices, the canonical 



form of matrices w and z is found as 

(w-i \ 
w\ 

w 2 

W 2 



-zi\ 
zi 



-Z2 

z 2 



(42) 



/ 



V ■ / V • 

where Wi = y/1 — n i} Zi = ^fnl. 

We now show that we can obtain the idempotent R 2 — R by applying 
the secondary Bogolyubov transformation (40), which automatically leads 
to Eq. (37) . Indeed, using the inverse transformation of (40) with matrices 
w and z given in Eq. (42), we found that the modified quasiparticle-density 
matrix Q can be obtained as 



W t Q W = 



where 



and 



W 



{y/T=n)< (VI) 



n) 1 — n 



It 



Q 



in 



WW t 



Qo 



«ata)) (M) 
«atat)) 1 - ((at a)) 




1 



Qo — Qo 



(43) 



(44) 



(45) 



due to Eq. (39). This result shows another way of deriving the modified 
quasiparticle-density matrix Q (34) from the density matrix Qo of the 
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modified quasiparticles (at, aj). This matrix Q is identical to the zero- 
temperature quasiparticle-density matrix Q (27). Substituting this result 
into the right-hand side (rhs) of Eq. (33), we obtain 

WQ U , 



where 



R 



u = wu 



(46) 



u* V* 
V u 



(VT^n)*u* + {Vn)*v 

~nV + VEU* 



(47) 



VT - nV + Vnu* VT 

This equation is the generalized form of the modified Bogolyubov coeffi- 
cients Uj and vj given in Eq. (38) of Ref. 9 From Eqs. (18), (44), and (47), 
it follows that UW = 1, i.e. transformation (46) is unitary. Therefore, from 
the idempotent (45) it follows that R 2 — R. 

Applying the Wick's theorem for the ensemble average, one obtains the 
expressions for the modified total energy £ 



where 



kl 



2^2 V ijklTkl 
kl 



(48) 



(49) 



From Eq. (46) we obtain the modified single-particle density matrix p and 
modified particle-pairing tensor f in the following form 



p = U T nU* + V\l - n)V + U 1 



\/n(l - n) 



V + V^y/nil ~n)U* , (50) 



f = U T nV* + V f (l - n)U + U 1 



\/n(l - n) 



U + V^n{l-h)V* . (51) 



As compared to Eq. (21) within the FT-HFB approximation, Eqs. (50) 
and (51) contain the last two terms <~ [yjn{l — and ~ \/n(l — n), 
which arise due to quasiparticle-number fluctuation. Also the quasiparticlc 
occupation number is now n [See Eq. (36)] instead of n (14). 

We derive the MHFB equations following the same variational proce- 
dure, which was used to derive the FT-HFB equations in Ref. 3 According 
it, we minimize the grand potential 6Q = by varying U, V, and n, where 



il = £-TS~fiN 



(52) 
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The MHFB equations formally look like the FT-HFB ones, namely (22) 

-1. J.) (2) -*(";)■ 

where, however 

H = T - p, + f (54) 

with f and A given by Eq. (49). The equation for particle number N within 
the MHFB theory is 

N = Trp . (55) 

By solving Eq. (53), one obtains the modified quasiparticle energy Ei, which 
is different from Ei in Eqs (22) due to the change of the HF and pairing 
potentials. Hence, the MHFB quasiparticle Hamiltonian £Tmhfb can be 
written as 

H - p,N W H M HFB = E Q + J2 ' ( 56 ) 

i 

instead of (2). This implies that the approximated density operator D (32) 
within the MHFB theory can be represented in the form similar to (12), 
namely 

V = Pmhfb = YllnJ^i + (1 - ni){l - M)]. (57) 

i 

From here it follows that the formal expression for the modified entropy S 
is the same as that given in Eq. (25), i.e. 

S = - X^ iln ^ + - ^)ln(l - Si)] , (58) 

i 

Using the thermodynamic definition of temperature in terms of entropy 
1/T = SS/SS and carrying out the variation over 5n i} we find 

— =Ei = T—=Tln[ — 

orii orii \ rii 

Inverting Eq. (59), we obtain 



(59) 



* - ^TT ' (60) 

This result shows that the functional dependence of quasiparticle occupa- 
tion number n 4 on quasiparticle energy and temperature within the MHFB 
theory is also given by the Fcrmi-Dirac distribution of nonintcracting quasi- 
particles but with the modified energies Ei defined by the MHFB equations 
(53). Therefore we will omit the bar over fii and use the same Eq. (14) with 
Ei replaced with Ei for the MHFB equations. 
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2.3. MBCS theory 

In the limit with equal pairing matrix elements Gij = G, neglecting the 
contribution of G to the HF potential so that f = 0, the HF Hamiltonian 
becomes 

"Hij = (e» - ■ (61) 
The pairing potential (49) takes now the simple form 

fc>0 

The Bogolyubov transformation (4) for spherical nuclei reduces to 

(-) j+m aj. m = Uj (-y +m aj - m v 3 a] m , (63) 
while the secondary Bogolyubov transformation (40) becomes 9 

(-Y +m a^ m = ^T^(-y+m a ._ m + ^ a t m (64) 

The U, V, 1 — n, n, and y/n(l — n) matrices arc now block diagonal in 
each two-dimensional subspace spanned by the quasiparticle state \ j) and 
its time- reversal partner \j) = (— Y +m \j — m) 



. 1 — nj \ ( rij 



x/ nf 1 - n) = T -vS-q-n,) 



(66) 



Substituting these matrices into the rhs of Eqs. (50) and (51), we find 

p. ~. = (1 - 2n 3 >? + nj - 2Jn j (l-n j )u j v j , (67) 



t. = -(1 - 2nj)«j«j + J n 3 -(l - nj)(«? - «*) . (68) 
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Substituting now Eqs. (68) and (67) into the rhs of Eqs. (62) and (55), 
respectively, we obtain the MBCS equations for spherical nuclei in the fol- 
lowing form: 

A = G fij [(1 - 2n>^- - yjn^l-nj)^ - v>)} , (69) 

3 

N = 2J2^ [(1 - 2n>f + nj - 2yJ nj (l -n^UjVj] . (70) 

3 

Comparing the conventional FT-BCS equations, we see that the MBCS 
equations explicitly include the effect of quasiparticle-number fluctua- 
tion ~ SjVj in the last terms at their rhs, which are the thermal gap 
— G^2j fljy/rij(l — rij)(uj — Vj), and the thermal- fluctuation of particle 
number SN = J2j &Nj = QjUjVj(6Afj) in Eq. (70). These terms arc 

ignored within the FT-BCS theory. Hence Eqs. (69) and (70) show for the 
first time how the effect of statistical fluctuations is included in the MBCS 
(MHFB) theory at finite temperature on a microscopic ground. So far this 
effect was treated only within the framework of the macroscopic Landau 
theory of phase transition. 5 



3. Phonon-damping model in quasiparticle representation 

The quasiparticle representation of the PDM Hamiltonian 13 is obtained by 
adding the superfluid pairing interaction and expressing the particle (p) and 
hole (h) creation and destruction operators, a\ and a s (s — p,h), in terms 
of the quasiparticle operators, a\ and a s , using the Bogolyubov's canonical 
transformation. As a result, the PDM Hamiltonian for the description of 
EA excitations can be written in spherical basis as 

H = E E 3 a ]m a jm + E U^XnM^ + 
jm Xfj,i 



Xfxi 



33' 



'33' 



l 33' 



A] f (X^+A n ,(\fi) 



(71) 



33 



Bjj> (^m) + Bjji (A/i) 




where A = v2A-j-T. The first term at the rhs of Hamiltonian (71) corre- 
sponds to the independent-quasiparticlc field. The second term stands for 
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the phonon field described by phonon operators, b\- and bx^i, with multi- 
polarity A, which generate the harmonic collective vibrations such as GDR. 
Phonons arc ideal bosons within the PDM, i.e. they have no fermion struc- 
ture. The last term is the coupling between quasiparticle and phonon fields, 
which is responsible for the microscopic damping of collective excitations. 
In Eq. (71) the following standard notations are used 

A h> ( A ^) = E (jmj'm'\Xfx)a] m a], m , , (72) 

mm/ 

B )r (V) =-£ (-y'- m ' (jmf m'\\n)a) m a jlm , , (73) 

mm' 

with < > (— ) A ~ M (A — /i). Functions u^} = UjVj> + VjUy and v"--} = 

UjUji — VjVji are combinations of Bogolyubov's u and v coefficients. The 
quasiparticle energy Ej is calculated from the single-particle energy ej as 



E j = ^/(e>-e F )Z + A* , e'j = ej - Gv^ , (74) 

where the pairing gap A and the Fermi energy ep are defined as solutions 
of the BCS equations. At T ^ the thermal pairing gap A(T) (or A(T)) 
is defined from the finite-temperature BCS (or MBCS) equations. 

The equation for the propagation of the GDR phonon, which is damped 
due to coupling to the quasiparticle field, is derived making use of the 
double-time Green's function method (introduced by Bogolyubov and 
Tyablikov, and developed further by Zubarev 14 ). Following the standard 
procedure of deriving the equation for the double-time retarded Green's 
function with respect to the Hamiltonian (71), one obtains a closed set of 
equations for the Green's functions for phonon and quasiparticle propaga- 
tors. Making the Fourier transform into the energy plane E, and expressing 
all the Green functions in the set in terms of the one-phonon propagation 
Green function, we obtain the equation for the latter, G\i(E), in the form 

where the explicit form of the polarization operator P\i{E) is 



33 



( V 3j' ] ) 2 ( n 3- n 3')( e 3- e 3') 



£2_ (e . +e .,)2 

(76) 
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The polarization operator (76) appears due to ph - phonon coupling in the 
last term of the rhs of Hamiltonian (71). The phonon damping "f\i(u)) (u) 
real) is obtained as the imaginary part of the analytic continuation of the 
polarization operator P\i(E) into the complex energy plane E = u)±is. Its 
final form is 



7«M = rrrEl/. 



33 



(uf 3 )Y(l - n 3 - ny )[6(E- E, - Ey ) - S(E + E 3 + Ey )] - (77) 



(«y ) 2 K" - n r )[5(E - E 3 - + Ey) - S(E + ej - ey)]j. 

The energy Q of giant resonance (damped collective phonon) is found as 
the solution of the equation: uj — w\i — P\i(Co) = . The width T\ of giant 
resonance is calculated as twice of the damping j\(uj) at u — Q, where A = 
1 corresponds to the GDR width Tgdr- The latter has the form 

Tgdr = YMvVf^ - U P- n h)S(E GDR -E p - E h ) + 



F 2 bL'^K' - n s)S(E GDR -E s + E s ,)\ 



(78) 



where (ss r ) = (pp r ) and (hh r ) with p and h denoting the orbital angular 
momenta j p and jh for particles and holes, respectively. The first sum at the 
rhs of Eq. (78) is the quantal width Tq, which comes from the couplings of 
quasiparticle pairs [aJ^a^LAf to the GDR. At zero pairing they correspond 
to the couplings of ph pairs, [a P ®aj^LM to the GDR. The second sum comes 
from the coupling of [a\ ® ol~\lm to the GDR, and is called the thermal 
width Tt as it appears only at T 7^ 0. At zero pairing they are pp (hh) 
pairs, [a J (g> ci~,}lm (The tilde ~ denotes the time-reversal operation). 

The line shape of the GDR is described by the strength function 
S'gdr.(w), which is derived from the spectral intensity in the standard way 
using the analytic continuation of the Green function (75) and by expand- 
ing the polarization operator (76) around to = -Egdr- The final form of 
SgdrM is 12 - 13 

c f ,,\ 1 7gdrH (7Q , 

t (w - £gdr) 2 + 7gdrM 
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The PDM is based on the following assumptions: 

al) The matrix elements for the coupling of GDR to non-collective ph 
configurations, which causes the quantal width Tq, are all equal to F\. 
Those for the coupling of GDR to pp (hh), which causes the thermal width 
r*r, are all equal to F2. 

a2) It is well established that the microscopic mechanism of the quantal 
(spreading) width Tq comes from quantal coupling of ph configurations to 
more complicated ones, such as 2p2h ones. The calculations performed in 
Refs. 15 within two independent microscopic models, where such couplings to 
2p2h configurations were explicitly included, have shown that Tq depends 
weakly on T. Therefore, in order to avoid complicate numerical calculations, 
which are not essential for the increase of Tgdr at T 7^ 0, such microscopic 
mechanism is not included within PDM, assuming that Tq at T = is 
known. The model parameters are then chosen so that the calculated Tq 
and -E/QDR reproduce the corresponding experimental values at T = 0. 

Within assumptions (al) and (a2) the model has only three T- 
independent parameters, which are the unperturbed phonon energy uj q , F±, 
and i*2- The parameters ui q and Fi are chosen so that after the p/i-GDR 
coupling is switched on, the calculated GDR energy -Egdr and width Tgdr 
reproduce the corresponding experimental values for GDR on the ground- 
state. At T 0, the coupling to pp and hh configurations is activated. The 
F2 parameter is then fixed at T = so that the GDR energy -Egdr does 
not change appreciably with T. 

4. Numerical results 

4.1. Temperature dependence of pairing gap 




012345 012345 



T (MeV) T (MeV) 

Fig. 1. Neutron pairing gap (a) and heat capacity (b) for 120 Sn as functions of T. Solid 
and dotted lines show the MBCS and BCS gaps, respectively. 
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Shown in Fig. 1 (a) is the temperature dependence of the neutron pairing 
gap A„ for 120 Sn, which is obtained from the MBCS equation (69) using 
the single-particle energies determined within the Woods-Saxon potential 
at T = 0. The pairing parameter G„ is chosen to be equal to 0.13 MeV, 
which yields A(T = 0) = A(0) ~ 1.4 MeV. Contrary to the BCS gap 
(dotted line), which collapses at T c ~ 0.79 MeV, the gap A (solid line) does 
not vanish, but decreases monotonously with increasing T at T > 1 MeV 
resulting in a long tail up to T ~ 5 MeV. This behavior is caused by the 
thermal fluctuation of quasiparticle number in the MBCS equations (69). 
As the result, the heat capacity [Fig. 1 (b)] has no divergence at T c , which 
is seen within the BCS theory. 

4.2. Temperature dependence of GDR width 

The GDR widths as a function of T for 120 Sn obtained within the PDM 
are compared in Fig. 2 (a) with the experimental data and the prediction 
by the thermal fluctuation model (TFM). 16 




1 2 3 4 5 5 10 15 20 25 30 



T (MeV) Ey ( MeV) 

Fig. 2. (b): GDR width r G DR as a function of T for 120 Sn. The thin and thick solid lines 
show the PDM results obtained neglecting pairing and including the renormalized gap 
A = [1 + 1/(5 N 2 ]A, respectively. The gap A includes the correction SN 2 = A(0) 2 + 
l/2)/[(ej — ef) 2 + A(0) 2 ] due to an approximate number projection. The prediction by 
the TFM is shown as the dotted line 16 ; (b): GDR strength function at T = 1.24 MeV. 
The dashed and solid lines show the results obtained without and including the gap A, 
while experimental results are shown as the shaded area. 

The TFM interprets the broadening of the GDR width via an adiabatic 
coupling of GDR to quadrupole deformations induced by thermal fluctu- 
ations. Even when thermal pairing is neglected the PDM prediction, (the 
thin solid line) is already better than that given by the TFM, including the 
region of high T where the width's saturation is reported. The increase of 
the total width with T is driven by the increase of the thermal width Tt, 
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which is caused by coupling to pp and hh configurations, since the quantal 
width Tq is found to decrease slightly with increasing T. 12 The inclusion 
of thermal pairing, which yields a sharper Fermi surface, compensates the 
smoothing of the Fermi surface with increasing T. This leads to a much 
weaker T-dependence of the GDR width at low T. As a result, the val- 
ues of the width predicted by the PDM in this region significantly drop 
(the thick solid line), recovering the data point at T = 1 MeV. The GDR 
strength function obtained including the MBCS gap is also closer to the 
experimental data than that obtained neglecting the thermal gap [Fig. 2 
(b)]- 



1.6 




1 2 3 0.4 0.8 1.2 1.6 

T (MeV) T (MeV) 



Fig. 3. (a): Pairing gaps for 120 Sn averaged over thermal shape fluctuations versus T. 
Lines with triangles and crosses are the usual BCS proton and neutron pairing gaps, re- 
spectively, while those with diamonds and squares denote the corresponding pairing gaps, 
which also include thermal fluctuations of pairing fields, (b): GDR widths for 120 Sn versus 
T. Open squares, triangles, and diamonds denote the widths obtained without pairing, 
including BCS pairing, and thermally fluctuating pairing field from (a), respectively. 

The results discussed above have also been confirmed by our recent cal- 
culations within a macroscopic approach, which takes pairing fluctuations 
into account along with the thermal shape fluctuations. 17 Here the free 
energies are calculated using the Nilsson-Strutinsky method at T ^ 0, in- 
cluding thermal pairing correlations. The GDR is coupled to the nuclear 
shapes through a simple anisotropic harmonic oscillator model with a sepa- 
rable dipole-dipole interaction. The observables are averaged over the shape 
parameters and pairing gap. Our study reveals that the observed quenching 
of GDR width at low T in 120 Sn and 148 Au can be understood in terms of 
simple shape effects caused by pairing correlations. Fluctuations in pairing 
field lead to a slowly vanishing pairing gap [Fig. 3 (a)], which influences the 
structural properties even at moderate T (~1 MeV). We found that the 
low-T structure and hence the GDR width are quite sensitive to the change 
of the pairing field [Fig. 3 (b)]. 
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5. Conclusions 

It has been shown in the present lecture that the MHFB and MBCS theories 
are microscopic approaches, which take into account thermal fluctuations 
of quasiparticle number. These large thermal fluctuations smooth out the 
sharp SN phase transition in finite nuclei. As a result, the thermal pairing 
gap does not collapse, but decreases monotonously with increasing temper- 
ature T, remaining finite even at T as high as 4 - 5 MeV. This non- vanishing 
thermal pairing gap keeps the width of GDR remain almost constant at low 
T (< 1 MeV for 120 Sn) when it is included in the PDM. In this way the 
PDM becomes a semi-microscopic model that is able to describe the tem- 
perature dependence of the GDR width in a consistent way within a large 
temperature interval starting from very low T, where the GDR width is 
nearly T- independent, to the region when the width increases with T (1 
< T < 3 - 4 MeV), and up to the region of high T (T > 4 - 5 MeV), where 
the width seems to saturate in tin isotopes. 
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